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ABSTRACT 



P^ ■ We present Spitzer Infrared Spectrograph (IRS) observations of the ultra- 

^ ■ luminous X-ray source (ULX) NGC 6946 X-1 and its associated nebula MF 16. 

This ULX has very similar properties to the famous Holmberg II ULX, the first 
ULX to show a prominent infrared [O IV] emission line comparable to those 

'ES ■ found in AGN. This paper attempts to constrain the ULX Spectral Energy Dis- 

tribution (SED) given the optical/UV photometric fluxes and high-resolution 
X-ray observations. Specifically, Chandra X-ray data and published Hubble op- 
^ ■ tical/UV data are extrapolated to produce a model for the full optical to X-ray 

f^^ ' SED. The photoionization modeling of the IR lines and ratios is then used to 

rn 



test different accretion spectral models. While either an irradiated disk model or 
an 0-supergiant plus accretion disk model fit the data very well, we prefer the 
O ' latter because it fits the nebular parameters slightly better. In this second case 

the accretion disk alone dominates the extreme- UV and X-ray emission, while an 



0-supergiant is responsible for most of the far-UV emission. 



H ■ Subject headings: black hole physics — galaxies: individual (NGC 6946) — in- 

■ ■ frared: ISM — X-rays: binaries 



1. INTRODUCTION 

The suggestion that ultraluminous X-r ay sources (ULXs) in nearby galaxies may harbor 



intermediate-mass black holes (IMBH, e.g. IColbert fc Mushotzkylll999l ) has major implica 



tions on the studies of black hole formation and evolution. However, the evidence for this 
new type of black hole has been slow to accumulate. Other explanations have been proposed 
to account for the high X-ray luminosities of ULXs, such as super-Eddington accretion 
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(iDone fc Kubota 2006 



2008 



Abolniasovll201ll ) or beaming (JKording et al 



: kuncic et al.l 120071 iKineJboOSi: ISoria fc KuncidboOSl : Berghea et al. 



2002 



Kin 



20081 ) ■ Som e galactic binaries 



are k nown to show sup er-Eddington luminositi es, such as GRS 1915+105 ( iFender fc Belloni 
20041 ) and V4641 Sgr f lRevnivtsev et al.ll2002l ). However conclusive evidence on the true 
nature of some or all known extra-galactic ULXs has remained elusive. 

One of the most interesting developments in ULX history is the discov ery in recent years 
of lar ge ionized bubble nebulae around some of the most famous ULXs (jPakull fc Mirioni 
20021 ). These are much larger and more luminous than regular supernova remnants. In 
addition, high-ionization lines have been de tected, such as He II A4 686 (LP. = 54.4 eV; e.g. 
Abolmasovet al.lboOTt IMoou et al.lboilh . IXaaret fc Corbel! J2009I ) detected [Ne V] A3426 
from the nebula associated with the ULX in NGC 5408 (LP. = 97.1 eV). The detection of 
these lines suggests photoionization of the surrounding gas by the hard radiation from the 
accreting ULX, similar to what happens with active galactic nuclei (AGN). 

It is also remarkable that there may exi st such an object in our own g alaxy. The first 
recognized microquasar, the famous SS 433 (IStephenson fc SanduleakI 119771 ) is possibly the 
"closest" relati ve to the nebulae- surrounded ULXs. It is surrounded by a 100 pc radi o 



nebula, W 50 ( iDubner et al.lll998l ) and is most likely super-Eddington (e.g. iFabrikal 120041 ). 



Indeed, SS 433 could be an example of both beaming and super-Eddingto n accretion, the 



combination of which could easily explain even the most luminous ULXs (JBegelman et al. 



20061 : iPoutanen et al.l 120071 ) . It is seen almost edge-on and therefore most of its strong U V 
and X-ray radiation is likely not visible. However, if such objects exist in nearby galaxies, 
those seen at small inclination angles should appear as bright ULXs. Recently an object very 



similar to SS 433 was discovered in the galaxy NGC 7793 ( jPakuU et al.l l2010l : ISoria et al. 
20101), showing radio lobes and hot spots, like a bona-fide quasar. 



Vazquez et al.l (120071 ) conducted the first infrared spectroscopy study of ULXs. They 



obtained high-resolution Spitzer spectra, of six ULXs in NGC 4485/90. IVazquez et al.l (120071 ) 
looked for the mid-infrared [O IV] high ionization line, which has an equivalent ionization 
potential to He II. However, apart from a possible weak detection in one object in their 
sample (ULX-1, which is probably a supernova remnant), [O IV] was not found in the 
NGC 4485/90 ULXs. 

The first reliable detect i on of [0 IV] associated with a ULX was the famous Holm- 
berg II ULX ( Berghea et al.ll2010a| ]bl). It is also one of the most interesting examples of 



an ionized nebula associated with a ULX (iPakuU &: Mirionil l2002l : iLehmann et al. 



2005 



Abolmasov et al.l 120071 ). He II is detected in a sma. l l nebular region ar ound the ULX, the 
"Heel" of the "Foot Nebula" JPakuU fc Mirionil |2002[ iKaaret et aJl2004[ ). An OB star is the 
likely counterpart. The morphology of He II emission and other optical lines is consistent 
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with the photoionization of the nebula by the ULX (jPakuU &: Mirionil 120021 ) . lKa.aret et al. 



Miller et aLrj2005[ ) detected 



( 120041 ) confirmed this using the Hubble Space Telescope (HST). 
extended radio emission from a ~50 pc diameter region coincident with the position of the 
ULX. These authors conclude that the flux and morphology of the radio emission is incon- 
sistent with emission from a SNR or an H II reg i on, an d is probably associated with the 
ULX activity. Using new observations, ICseh et al.l ( 12012| ) suggest that the radio emission is 
probably caused by a weakly collimated outflow, rather than a jet. 

The ULX which we investigate in this paper (which w"e wil l call NGC 6946 X-1 from 
now on) was first identified as such by [Roberts fc ColbertI (120031 ) . who showed that the X- 
ray variability of this source indicates a binary origin. MF 16 has been known as a very 
lumino us SNR iMatonick fc FesenI (119971) associated w ith a ROSAT X-ray source and a radio 
source (Ivan Dyk et al.lll994J ). lAbolmasov et al.l ( 120081 ) used optical spectroscopy to conclude 
that the ULX is most likely responsible for the high ionization lines detected in the MF 16 
nebula. They suggest that the ULX is probably bright in UV/EUV and not only in X-rays 
(p redicted absorbed fl ux at 1000 A: 10~^^ erg cm~^ s~^ A ~^). This was actually verified 
by lKaaret et al.l ( l2010l ). who used the ACS camera on the Hubble Space Telescope to detect 
far-UV emission from the ULX for the first time. The data does not allow the authors to 
conclude with certainty what the UV source is, but they suggest it is a combination of disk 
emission and contribution from an OB supergiant. 



Recently a tentative quasi-periodic oscillation of ~8.5 MHz was detected (JRao et al. 



20101) using XMM-Newton X-ray data which - if confirmed - makes the NGC 6946 X-1 
the fourth ULX with such timing features (2 ULXs in M82 and NGC 5408 X-1) NGC 
6946 X-1 is very similar to the Holmberg II ULX in many respects. iBerghea et al.l (l2010bl ) 
used Mappings HI photoionization and shock modeling to show that NGC 6946 X-1 and 
Holmberg II ULX high-ionization line fluxes primarily result from photoionization rather 
than shocks (see also Section 4.1.). 

In this paper we use IR Spitzer spectroscopic observations to constrain the UV and 
soft X-ray emission of the NGC 6946 X-1. The IR observations, when combined with 
previously published optical, UV and X-ray data, suggest that the powerful UV and soft 
X-ray emission from NGC 6946 X-1 photoionizes the surrounding medium. In Section 2 we 
present the Spitzer IRS data analysis and Chandra X-ray spectral fitting results. In Section 3 
we use our data and other published data to construct an SED of the ULX which is then used 
as photoionization source in CLOUDY modeling. Finally, in Section 4 we use our results to 
compare the ionization properties of this famous ULX with those of standard AGN. 
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DATA PROCESSING AND RESULTS 



2.1. Spitzer mS Observations 



The data presented here are archival Spitzer data, taken on 2009 February 21. The IRS 
spectral map is described by four parallel pointings and 12 perpendicular at half-slit steps, 
for the Short High (SH, in the range 9.9-19.6 /im) module. The Long High map (LH, in the 
range 18.7 — 37.2 //m) has three parallel pointings and 7 perpendicular. 

The total exposure time was 1.4 ks and 1.3 ks for SH and LH observations, respectively. 
The data were preprocessed by the IRS pipeline (version 18) at the Spitzer Science Center 
prior to download. Spectral maps were constructed using "BCD-level" proce ssed data in 
coniunction with CUBISM (CUbe Builder for IRS Spectra Maps) version l.E^J ( ISmith et al. 



2007; Kennicutt et al. 



i s ^25 % fISmith et al 



20031). The absolute calibration error on fluxes obtained with CUBISM 



200. 



^. We adopt a distance of 5.1 Mpc to NGC 6946 as in lKaaret et al. 



( 120101 ) (1^^ is 25 pc). The LH maps obtained with CUBISM have intrinsic sizes of 26.4" by 
17.6" (660 pc by 440 pc), while the SH maps are 23" by 13.8" (575 pc by 345 pc). However, 
in order to avoid aperture effects, the spectra presented here were extracted from a like-sized 
aperture. Thus the SH and LH spectral data result from an equally sized square aperture of 
size 13.3" (332.5 pc), which is composed of 9 LH pixels centered on the ULX position. 

The high-resolution spectra (10 /xm— 37 /xm) are plotted in Figured] Since the spectra 
are extracted from matched apertures, the spectral lines in both the SH and LH modules 
are properly calibrated and the continuum is properly matched. We used SMART (Spec- 
troscopic Modeling Analysis and Reduction Tool, see iHigdon et al.l |2004[ ) version 8.1.2 to 
measure line fluxes. The fluxes of detected lines and upper limits for non-detections are 
prese nted in Table [H We no te that the [O IV] line is detected at a signal-to-noise ratio 
> 12. lAbolmasov et al.l ( 120081 ) estimates an extinction of Ay = 1 . 54 ma g in the direction of 



NGC 6946 X-1. For the infrared lines, we estimate from iDraind (Il989[ ) that the extinction 
is Aa < 0.1, and therefore negligible. 

The SH and LH maps for the detected lines are presented in Figure [2l They are shown 
together with two Hubble Space Telescope (HST) images of the MF 16 nebula in UV and 
Ha emission. A nearby (200 pc away) star-forming region and the associated HII region 
can be seen just outside the Spitzer aperture. As can be seen in the Spitzer spectral maps, 
the low-ionization lines ([Ne II], [S HI], and [Si II]) are contaminated by the surrounding 
star-forming region. In contrast, the higher ionization lines [O IV] and [Ne HI] are clearly 



^Available at: http://ssc.spitzer.caltech.edu /archanaly/contributed/cubism/ 



-5- 



concentrated at the NGC 6946 X-1 location, suggesting a source of high-energy photons. 



2.2. Archival Data 



2.2.1. X-ray spectrum: 

We use X-ray spectroscopy in conjunction with published optical/UV data to model the 
X-ray emission from the ULX and to extract the full optical to X-ray spectrum. We also use 
it to estimate the UV emission by extrapolating the model to lower energies. To model the X- 
ray spectrum we used the best available data from the longest (~58 ks) C/iandra observation 
in the archive (Obs. ID 1043, taken on 2001 September 7). The data were processed by 



Berghea et al.l (120081 ) (U45 in this paper) and the spectrum is plotted in Figure O 



We note that iFridriksson et al.l (120081 ) find long-term X-ray flux variability of at most 
30% in NGC 6946 X-1 based on Chandra, XMM-Newton and ROSAT data. However for 
completeness we have chosen the best data available to represent the source (e.g. Chandra 
data rather than ROSAT data), as well as one in which the luminosity is approximately 
average over the 10 years of observations (i.e. the 2001 data set). This average X-ray 
spectrum and luminosity best represents the ULX in a state in which the IR-line fluxes can 
be modeled. In addition, the equilibrium time-based recombination for the high ionization 
lines are approximately 3000 to 100 years for He II and [O IV] r espectively— well above the 
variability time-scales of the X-ray source (JBerghea et al.ll2010bl ). 



2.2.2. HST UV and optical photometry: 



We use the measurements obtained by iKaaret et al.l (120101 ) to constrain the optical and 
UV emission from the accreting source and its companion. The B, V, I and UV (F140LP 
filter) extinction-corrected fluxes are plotted as circles in Figure [31 



CLOUDY MODELING 



3.1. Fitting the X-ray Spectrum and Optical Data 

We used XSPEC version 12.4 to model the spectrum between 0.3 and 8.0 keV. The 
X-ray spectrum is the most important component of the total SED since it is likely the 
only observed emission that directly emanates from the ULX, and is not contaminated by 
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the surrounding star-forming region. For this reason we have taken many precautions in 
modehng the X-ray spectrum, however we fit the X-ray and optical/UV hterature data 
simultaneously. In this section we describe the models we chose to fit the X-ray spectrum 
and extrapolate to the unknown UV region not covered by literature data. 

Probably the most common and simplest model for ULX X-ra y spectra is the two- 
component model with a multi-color accretion disk component (MCD. lMitsuda. et al.lll984j ) 
and a power-law for the hard tail. However, this model c annot be extrapolated to UV 
where the source of IR-lines originates (JDone fc Kubotall2006l . e.g.), and therefore we use the 
physically motivated Comptonization models described in detail below. 



For all models, following iBerghea et al.l ( 120081 ). the Galactic absorption column density 
was fixed at 2.02x10^^ cm~^, and an additional local abso rption (Nh) model was ad ded. 
We assume that the metallicity is approximately Solar (see iGarcia-Benito et al.l |2010| . and 
references therein). 

DISKIR Model: The DISKIR model in XSPEC is a particularly useful Comptoniza- 
tion model to f it both optical and X -ray data, because it also includes emission from an 
irradiated disk (JGierlihski et al.ll2009[ ). We found that the best fit model has a disk of radius 
500 Rjn (the inner disk radius), and the irradiated flux is at 1.5% of the total flux. The 
fraction of the flux thermalized in the inner disk was fixed at 0.1. In addition, the electron 
temperature kTe and the radius of the Compton illuminated disk r,> T- were set at 100 keV 
and 1 .1, respectively. These values are typical for this model (e.g. 
20081 : iGrise et al.ll2012h . 



Gierlihski et al. 2009 



BMC Model: We also used th e Bulk Motion Compt onization model (BMC in XSPEC), 
described by lTitarchuk et al.l ( 1l997l ). iGliozzi et al.l (120111 ) have suggested that this model is 
useful for ULX-like black-hole accretion scenarios. The best fit is obtained for a fraction of 
Comptonized flux of 0.12%. 

SIRF Model: Finally we tried the Self-IrRadiated F unnel model (SIRF in XSPEC), 
which describes the emission from a "supecritical funnel" (jAbolmasov et al.ll2009l). This is 
create d in objects similar to SS 433, by the wind formed at high accretion rates. lAbolmasov et al, 
koO^ obtained a reasonable fit for XMM-Newton data of NGC 6946 X-1. The funnel open- 
ing angle was fixed at 20 degrees, the inclination angle was set to zero, and the accretion rate 
fixed at an Eddington ratio of 10^. These model parameters were set based on the best-fit 
results to the X-ray dat a. We obtained similar results for inner temperature and radius as 
Abolmasov et al.l ( 120091 ) . but unfortunately could not find an acceptable fit using this model 
as Table |2] shows. As a result, this model has been disregarded as a viable model to fit the 
NGC 6946 X-1 SED in the discussion to follow. 



-7- 



The results from the X-ray fits are presented in Table [2l and the models are plotted in 
Figure [31 The first two models, DISKIR and BMC both provide equally good fits to the X- 
ray data and will be used to construct the input Spectral Energy Distribution (SED) for the 
CLOUDY modeling. However, as seen in Figure [31 they predict very different UV and optical 
emission, and therefore should yield very different results in the photoionization simulations. 
The physical difference between these two models is that the BMC model requires a bright 
optical companion to fit the optical and UV data. As seen in Figure [3l the DISKIR model 
fits the optical and UV data equally well with a bright (supergiant) companion (dotted red 
line) or an irradiated disk with a less luminous companion (solid red line) . 

The r nodel for the 09 I su pergiant shown as the solid line in blue in Figure [31 is a 
TLUSTY (JHubeny &: Lanzlll995l ) O star model with effective temperature 31000 K, log g = 3.5, 



and a luminosity of 3.16 x 10 erg s . This fits the optical HST data (round black circles 
in Figure [3]) very well. We add this star to both DISKIR and BMC models (dotted lines), 
but we remind our reader that the irradiated disk alone (solid red line) fits the optical-UV 
data equally well. Finally, we find that the SIRE model peaks just below the [O IV] edge, 
and is faint in the optical-UV, and therefore does not provide a good fit for this data as 
described above. 



3.2. CLOUDY Photoionization Modeling 

We use DISKIR and BMC SED models shown in Figure [31 as input spectra for CLOUDY. 
The bolometric luminosities that result from these models are 8.02 and 2.51 x 10^*^ erg s~^ 



respectively. We assumed a Solar r netallicity as per iGarcia-Benito et al.l (I2OIOI). Previous 



optic al studies of NGC 6946 X-1 (JDunne et al. 



200' 



I 



Blair et al. 



2001 



Abolmasov et al. 



20081 ) indicate that the ionized gas from the nebula was swept inside a thin shell 
width) around the ULX. 



1 pc in 



Estimates based on line ratios of [SII] suggest that the derisity o f the gas in the shell is 
quite high (400—600 cm~^) (iBlair et al.ll200ll: [Abolmasov et al.ll2008[) . If we use the ratio of 
the [S III] IR lines to estimate the density (iDudik et al.ll2007|). we obtain densities lower than 
150 cm~^ for temperatures in the range 10,000—20,000 K ([Abolmasov et al.|[2008l ). However, 
Figure [21 shows that the [S HI] lines are contaminated by the nearby star- forming region and 
therefore we prefer the higher value obtained from previous optical line ratios estimates, 
since these were obtained from higher resolution maps. T hus we assume a spherical shell 
geometry, a density of 400 cm~^, and a filling factor of 0.2 ( IBlair et al.|[200l[ ). 



The inner and outer radii of the shell were varied to obtain the best fit for the higher 



ionization lines [O IV] and [Ne III], since these IR-lines are less contaminated by surrounding 
star-forming regions as illustrated in Figure [2] We also u se the previously measured He II 



line flux (1.3 x 10 erg s ) as a constraint (IWong||2010l ). since this line has an ionization 



potential similar to [O IV] and [Ne III] and likely emanates from a similar region to these 
IR-lines. The predicted lines from the CLOUDY simulations are shown in Table |3] and 
Figure HI 

Previous optical observations have s hown that radius of the shell (or shells) is 7—10 pc 
and the thickness is of the order of ~1 pc ( iBlair et al.ll200ll ). The DIS KIR model i s cons istent 
with a shell of radius 20 pc and is thin compared to the findings of iBlair et al.l (l200l[ ) by a 
factor of ^1 0. The BMC mo del predicts radial dimensions that are closer to those values 
observed bv lBlair et al.l (2001), bu t the shell predicted is smaller by a factor of ~1.5— 2. A 
larger shell than iBlair et al.l (|200l[ ) is expected for the DISKIR model, because it provides 
more ionizing photons. Both models are very good fits to the high ionization IR lines as 
shown in Figure HI However, the BMC model over-predicts most low-ionization lines, while 
DISKIR under-predicts them. 

We finally note that both models under-predict the amount of gas compared to the 
Hydrogen column densities obtained from spectral fits. This could be explained if part of 
the absorption takes place very close to the accretion disk and not near the line-emitting 
region. 



DISCUSSION 



4.1. The origin of the [O IV] emission 



Berghea et al.l ( l2010bl ) discuss in detail possible scenarios for the [O IV] detected in the 
Homlberg II ULX in relation with AGN and star-forming regions. The conclusions of that 
paper apply to NGC 6946 X-1 as well. The [O IV] is clearly associated with the ULX and not 
with the star-formation nearby (Figure [2]). Moreover, the hottest stars alone cannot produce 
th e [O IV] measu r ed from the MF 16 nebula. For more details we again direct our reader 



to iBerghea et al.l (|2010bl ) and the references therein. An interesting result of that paper 
was that hot stars can have a ir nportant effect zf a hig h energy source is also present. The 
detailed CLOUDY simulation in iBerghea et al.l (l2010bl ) showed that an O star can increase 
the [O IV] luminosity by a factor of 3. This is also true in our case if the ULX companion 
is an 0-star. We draw a similar conclusion based on the He II line, which has a very similar 
ionization potential to [O IV], and therefore requires a very high energy source. 



The results of our CLOUDY modeling cannot easily discriminate between an 0-star 
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plus accretion disk and an irradiated disk (Figure [3]). As discussed in Section 3.2, previous 
opt ical data suggest that the radius of the shell is 7-10 pc and the thickness is on order 1 
pc ( iBlair et al.ll200l[ ). Table [3] shows that the BMC model predicts a radius and thickness 
that is much closer to the optical findings, than the DISKIR predictions. Finally, Figure H] 
show that the BMC model predictions for He II and [Ne III] are very slightly better than the 
DISKIR model predictions. For all of these reasons, we prefer the BMC model (the 0-star 
plus accretion disk model) to the DISKIR model (an irradiated disk) to describe NGC 6946 
X-1. 



Berghea et al.l (l2010bl ) explored shocks as a possible mechanism for the lines for both 



NGC 6946 X-1 and Holmberg II ULXs. They found that the He II line is probably produced 
by shocks in the case of another famous nebula-surrounded ULX in Holmberg IX. However for 
NGC 6946 X-1 and Holmberg II ULXs the He II emission is consistent with photoionization 
and not with shocks. To further confirm that the [O IV ] alone is not dominated by shocks, 
we reproduced Figure 10a from iBerghea et al.l (l2010bl ) for NGC 6946 X-1. We used the 
brightest pixel in the [O IV] spectral map to estimate a peak emission for this line. We 
note that this is in fact a lower limit for the estimate, since it represents the average fiux 
over the large LH pixel. As shown in Figure the shock v elocity requir e d to p roduce the 



peak [O IV] in NGC 6946 X-1 is in excess of 300 km s"^. iDunne et all ( 120001 ) find shock 
velocities from Ha are at most 270 km s~^. This suggests that the shocks needed to produce 
the high- ionization lines [O IV] and He II are simply not present, based on the Ha data. 

The ultimate source of X-ray and high-ionization emission in ULXs has been widely 
debated. Indeed, multiple mechanisms have been proposed to explain the high-energy emis- 
sion, including beamed stellar-mass black holes, super-Eddington stellar mass black holes, 
and finally, intermediate mass black holes (IMBHs). Evidence for one or more of these sce- 
narios has been found for a variety of ULXs. While our observations and modeling of the 
ULXs in Holmberg II and the M16 nebulae strongly suggest these objects are not super- 
nova remnants or beamed stellar-mass black holes, this analysis cannot distinguish between 
a super-Eddington, stellar-mass black hole or an IMBH as the source of the X-ray, opti- 
cal, and high ionization emission. This is due in part because the X-ray continuum is fit 
equally well with models describing either scenario. We note on the other hand, that the 
infamous SS 433 in our Galaxy, is a known beamed, super-Eddington, stellar-mass black 
hole with a clear host nebula. All evidence provided in the literature points to shocks or 
winds/outflows as the primary mechanism(s) for ionization in this object and the nebular 



regio n immediately surrounding the ULX (e.g. lAbolmasov et al.ll2010t iBegelman et al.lll980 



20061 ). Because SS 433 is the only super-Eddington stellar mass black hole with a nebula 
definitively known to date, this is the only such source we have to compare with NGC 6946 
X-1. Our analysis suggests that unlike SS 433, NGC 6946 X-1 has high ionization lines 
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(specifically the Oxygen line here and the He II from lAbolmasov et al.ll2008l . that are most 
likely not produced by shocks). In this singular way the famous SS 433 super-Eddington 
stellar mass black hole is very different from NGC 6946 X-1. 

The deficiency of super-Eddington sources in nebular environments precludes a detailed 
comparison of NGC 6946 X-1 with such objects. However, multiple active black holes in 
dense environments that produce high-ionization lines via photoionization are known to 
exist and are well documented in the literature. These sources, while not IMBHs, may prove 
an excellent comparative resource. For this reason, in the subsequent analysis, we focus 
attention on this comparison. 



4.2. Diagnostic line ratios 

The infrared [O IV] line (together with [Ne V]) is a signature of high ionizatio n and 
excitation usually associated with actively accreting black holes (e.g. iLutz et al.lll998l ). This 
line is often used to disentangle the starburst emission i n composite type galaxies, and is 



a good indicator of the AGN power (IGenze 



2004 



Smith et al.ll2004j : ISatyapal et al. 



2007 



et al 



1998 



Sturm et al. 



2002 



Satvapal et al. 



Melendez et al.ll2008r ). iGenzel et al.l fll998[ ) and 



Sturm et al.l (120021 ) have used the ratio of [O IV] /[Ne II] to distinguish between starbursts 
and AGN dominated galaxies in the mid-IR . In Figure E^ we c ompare the NGC 6946 X- 
1 mid-IR line ratios with data published by iGenzel et al.l ( 19981 ). The comparison clearly 
indic ates high excitation fo r the ULX, similar to those seen in AGN. The Holmberg II ULX 
from iBerghea et al.l (l2010al ) is also plotted here, though the [Ne II] flux for this ULX is an 
upper-limit. Both sources appear as AGN-like when this diagnostic is used. 

Other infrared line diagnostics have been proposed by iDale et al.l (120061 . 120081 ) to dis- 
entangle AGN-type ionizing sources from star formation, using the ratios [S HI] 33.48 Atm / 



[ Si II] 34.82 /im and [Ne HI] / [Ne II]. Figure |6]d shows the four regions defined by lDale et al. 



(J2006|) from their Figure 5, and their published data for different sources in the SINGS sam- 
ple. NGC 6946 X-1 is located well wit hin region II, together with most Seyfert-like nuclei. 
Holmberg II, is also plotted here from iBerghea et al.l (l2010al ). and shows ratios similar to 
low-luminosity nuclei, but we note that both [S HI] and [Ne II] are upper-limits for this 
object. We also note that the [S HI], [Si II], and [Ne II] are heavily contaminated by star 
formation as is clear in the spectral maps in Figure | 21 However , we also note that the 
AGN and star forming regions presented by lDale et al.l (|2006l . |2008|) are subject to the same 
contamination biases and aperture effects as the ULXs. We therefore still find it useful to 
compare the two ULXs with these sources. 
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The fact that the bale et all ( 120061 . boOSl ) ratios for NGC 6946 X-1 are Seyfert-hke and 
that the [O IV]/[Ne II] ratio is quite high, suggests that NGC 6946 X-1 is very gas rich 
compared to r nost ULXs and has a very strong UV source producing the ionizing radiation. 
We note that IVazquez et al.l (120071 ) found similar results in their Spitzer IRS study of the 
NGC 4485/90 ULXs. However they did not find clear evidence for [O IV] in their sample of 



object s. In examining this data further we find slightly discrepant ratios from I Vazquez et al. 
( 120071 ). From our analysis, only ULX-2 in their sample actually falls marginally inside 



region I, the other five ULXs being located in regions II and III. In communicating with the 
authors we find that the discrepancy is likely do to differences in reduction pipeline used for 
processing or the apertures used in either analysis. 

Finally we compared the mid-IR line ratios for NGC 6946 X-1 (^ar id Holmberg I I ULX ) 
with the those of other objects presented in the diagnostic diagrams of IWeaver et al.l ( 120101 ). 
The error bars for the ULXs in this plot are approximatel y the same size as t he data points 



and re present 25% absolute calibration error for CUBISM (jSmith et al.ll2004j ). IWeaver et al. 



(120101) find that high- ionizati on AGNs from the Swift Burst Alert Telescope (BAT) sam- 
ple ( jTueller et al.l l2008l 120101 ) distinguish themselves from other, less luminous AGN and 
starburst galaxies in their [Ne III]/[Ne II] and [O IV] /[Ne II] ratios. A replication of the 
Weaver et al.l ( 120101 ) plot is show in Figure [71 These high-ionization AGN have ratios very 
close to or over 1. In Figure [7] we over-plotted the NGC 6946 X-1 as well as the Holmberg II 
ULX — the only two ULXs known to produce prominent [O IV] emission. Interestingly, 
NGC 6946 X-1 occupies the overlapping region between the high ionization BAT AGN and 
their lower-luminosity counterparts. The results for NGC 6946 X-1 are consistent with pho- 
toionization which is similar to that seen in LINER or Seyfert-like nuclei rather than star 
forming in nature. 

All three diagnostic diagrams suggest that the [O IV] emission in Holmberg II ULX and 
in NGC 6946 X-1 is dominated by the ULX and that the surrounding star forming regions 
provide little contamination to this line. Bot h ULXs have been shown to be candi date in- 
termediate mass black holes by other authors (JLehmann et al.ll2005l : iRao et al.ll2010l ). While 
these findings cannot definitively prove that NGC 6946 X-1 is either an IMBH or a super- 
Eddington accretor, it does suggest that the environment may be very similar for this source 
and AGN-like objects. We finally mention that the correlation of NGC 6946 X-1 with ultra- 
luminous infrared galaxies (ULIRGs) in Figure [6^ and low ionization nucl ear emission line 



regions (LINERs) in Figure [TJ may lend some credence the predictions of iMcKernan et al. 



(120101 . I2OIII ) who sug gest that many LINERs might host ULX-like IMBHs. Ind eed many 
ULIRGs are LINERs JSatvapal et"aDl2007l : IPudik et al.ll2007l : ICarrillo et al1ll999h . 



-12- 



4.3. Black Hole Mass Estimates 

Whether NGC 6946 X-1 is a super-Eddington source or an IMBH cannot be definitely 
determined from the observed and modeled data presented in this paper. However as dis- 
cussed in Section 4.2, the ratios from NGC 6946 X-1 look very similar to the ratios of 
well-known AGN. We therefore find it useful to estimate the possible black hole mass of 
this ULX, based on AGN criteria with the caveat that the criteria assumes NGC 6946 X-1 
behaves much like an AGN which is not definitely known from the observed and modeled 
data presented here. The bullets below show four black hole mass estimates derived for 
NGC 6946 X-1. 



1. iGliozzi et al.l ( 1201 ll . |2009[ ) have proposed a new method to estimate BH masses using 
the BMC model. This method uses the photon index of the Comptonized part of the 
spectrum and the normalization to scale from well-known ste llar-mass BHs with well- 
define d spectral transitions. Using a scaling method defined by lShaposhnikov fc Titarchuk 
(120091 ) that is based on a detailed study of 8 differ ent Galactic black hole systems with 
17 spectral transitions, iGliozzi et al.l (120091 1201 ll ) extends these scaling techniques to 
extra-galactic black holes sources. Indeed one of these sources is a newly discovered 
ULX in NGC 3621. Thus this technique can be used in the case of NGC 6946 X-1 
as well, with the assumption that the states of ULXs mimic those of GBHB. It as- 
sumes that BHs of all masses show similar spectral states, at least with respect to 
Co mptonization . If we scale the results from our BMC fits to XTEJ1550- 564 (as done 
bv IGliozzi et al.ll201ll i n Figure 1) and use M© =9.1 & d=4.4 kpc from lOrosz et al. 
(|201l|), and 6 =72 from lciiozzi et al.l J201l[ ). we obtain a BH mass of «1,800 M©. 



2. Both DISKIR and BMC X-ray models predict a cool disk. Assuming a BH accreting at 
the Eddington limit, the mass of the BH can be estimated at 1 — 1.5x10^ M©. However, 
this estimate assumes that the ULX is in a typical high state. If this is not the case, 
the inner disk tempe rature cannot be used reliably to estimate the mass of the BH 
Jstobbart et al.ll2006l ). 



3. Another estimate can be made just based on the bolometric luminosity calculated 
from our models. Using the values listed in Table |2l and assuming again accretion at 
Eddington limit, we obtain masses of 130 and 450 M©. 



4. Finally, iDasyra et al.l (|2008[ ) found a good correlation between the luminosity of the 
[O IV] line and the mass of the central BH in standard AGN. We estimate a BH mass 
of 2.7x10^ M0 if this correlation is extrapolated to NGC 6946 X-1. 
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In conclusion these four methods predict a very wide range of black hole masses for 
NGC 6946 X-1. Indeed the estimates predict masses that span three orders of magni- 
tude: 10^ to IO^Mq, so emphasis on any one estimate is difficult to justify. However 
it is interesting to note that all four put this object within the IMBH regime for the 
given criteria (e.g. accretion at the Eddington rate). 



5. CONCLUSIONS 

Using the mid-IR S';)zt2;er spectral maps of NGC 6946 X-1 combined with Chandra X-ray 
and HST optical data from the literature, we model the spectral energy distribution of this 
unique ULX. We find the following major conclusions from our modeling and observational 
results: 



The SED constructed based on X-ray and optical spectral modeling indicates that 
NGC 6946 X-1 is well fit with two Comptonization models: The DISKIR model and 
the BMC model, both available in XSPEC. Spectral modeling did not permit a good 
distinction between an irradiated disk or a bright companion plus accretion disk as 
source of the optical emission. However, the nebula parameters (namely radius and 
thickness) predicted by C LOUDY using th e BMC model better match the findings 
from optical observations ( iBlair et al.l 120011 ) and the lines are slightly better matched 
with the BMC model in this case. We therefore prefer the BMC model to the DISKIR 
model. 

The detected [O IV] is consistent with photoionization by the strong X-ray emission 
from the ULX, and not with shocks or nearby star-formation. This suggests that the 
X-ray emission is not strongly beamed. The data cannot however discriminate between 
super-Eddington accretion and the IMBH scenario. 

We compared our results with known super-Eddington sources (i.e. SS 433) as well as 
standard AGN. We find that the emission line properties of NGC 6946 X-1 are not well 
matched to the SS 433 super-Eddington source, since the lines from this source have 
been shown to be produced primarily by shocks. We also compared NGC 6946 X-1 line 
ratios with multiple AGN diagnostics. The line ratios from NGC 6946 X-1 consistently 
show emission very simil ar to most low- luminosi t y AGN. We plott ed di agnostic ratios 
from 



Genzel et al.l (Il998f l and lSturm et al.l (l2002f l . iDale et al.l f l2006h . and 



Weaver et al. 



(2010|). All plots show this ULX falling safely within the AGN regime rather than star 
forming regions. This suggests that NGC 6946 X-1 may have photoionization and 
nebula properties that are similar to those seen in AGN. 
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• A variety of approximations were used to estimate the black hole mass from the ob- 
served properties of this ULX and resulted in a predicted black hole mass range span- 
ning three orders of magnitude. Assuming this ULX is accreting at the Eddington 
limit we find that NGC 6946 X-1 has a black hole mass between 10^ to 10^ Mq. 
Our calculations suggest that NGC 6946 X-1 could have a black hole mass within the 
IMBH regime, with the caveat that the criteria for estimating black hole mass assumes 
NGC 6946 X-1 behaves much like an AGN which is not definitely known from the 
observed and modeled data presented here. 



We are extremely grateful to Marcio Melendez for all of his help interpreting and plot- 
ting Figure [71 We are also very grateful to the referee whose comments have greatly improved 
the quality and clarity of this paper. This research has made use of the NASA/IPAC Ex- 
tragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National Aeronautics and Space Adminis- 
tration. This work is based on observations made with the Spitzer Space Telescope, which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology under a 
contract with NASA. SMART was developed by the IRS Team at Cornell University and is 
available through the Spitzer Science Center at Caltech. 
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Table 1. Measured infrared lines 



Line 


Flux 


S/N Ratio 


L 


(1) 


(2) 


(3) 


(4) 


[No II] 12.81 Mm 


2.5±0.3 


13.5 


7.9±0.9 


[Nc V] 14.32 /jm 


<0.5 




<1.5 


[Ne III] 15.56 Aim 


2.3±0.2 


18.9 


7.1±0.7 


[S III] 18.71 Mm 


1.3±0.2 


12.5 


4.1±0.7 


[Ne V] 24.32 ^m 


<2.1 




<6.7 


[O IV] 25.89 Mm 


1.9±0.2 


12.8 


6.0±0.5 


[Fe II] 25.99 Mm 


<1.2 




<3.7 


[S III] 33.48 Mm 


3.0±0.2 


24.3 


9.4±0.6 


[Si II] 34.82 Mm 


2.4±0.4 


9.90 


7.4±1.4 



Note. — Fluxes are in 10~ W cm" , luminosi- 
ties (L) in 10 erg s~ . If the measurements errors 
are smaller than the absolute calibration accuracy 
of 25%, the latter were used. For nondetcctions wc 
show 3ct upper limits. 



Table 2. X-ray model fits 



Model 


Nh 


kTi„ 


r/Ri„ 


Ax /dof 


L 


(1) 


(2) 
(10^^ cm~^) 


(3) 
(keV) 


(4) 


(5) 


(6) 
{lO^O erg 



DISKIR 

BMC 
SIRF 



4.26±0.11 0.118±0.003 
3.32+°-^^ 0.110±0.008 

-^■°^-o.os -'■*^-o.ii 



2.44±0.08 
1+0.08 
■ 0.10 



2.40^ 



1.16/183 
1.17/182 

1.85/182 



Note. — (1): X-ray model, as described in Section 3 .1. (2): Intrin s ic hy - 
drogen column density. The Galactic column density from lBerghea et al.| i 20081 ^ 
(2.02x10 cm~ ) was added separately. (3): Inner disk temperature. (4): Pho- 
ton index for DISKIR and BMC models a nd the inner radius in 1 0~ "spherisa- 
tion radius" units for the SIRF model (see |Abolmasov et al.|2009| ). (5): Reduced 
X values for the fit and the number of degrees of freedom. (6): Unabsorbcd (in- 
trinsic) luminosities between 0.1 and 10 kcV. 



Table 3. CLOUDY modeling results 



Model 


R,adius 


Depth 


He II 


[Ne II] 


[No V] 


[Ne III] 


[S III] 


[No V] 


[CIV] 


[Fe II] 


[S III] 


[Si II] 








A4686 


12.81 


14.32 


15.56 


18.71 


24.32 


25.89 


25.99 


33.48 


34.82 


(1) 


(2) 


(3) 
(pc) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


BMC 


5.1 


2.6 


37.14 


37.16 


36.97 


37.93 


38.02 


37.04 


37.82 


37.86 


38.24 


38.25 


DISKIR 


20.7 


0.06 


37.40 


35.98 


36.38 


37.69 


37.68 


36.45 


37.79 


36.57 


37.79 


37.29 



Note. — (1): SED model, as described in Section 3.1. (2): Inner shell radius. (3): The thickness of the shell. 
Line luminosities are in units of log erg s 
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Fig. 1. — Full IRS spectra obtained from matched apertures with CUBISM from the aperture 
shown in Fig. [2] 
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Fig. 2. — HST ACS images of the MF 16 nebula in NGC 6946 (left), and the Spitzer spectral 
maps of the infrared detected lines (right). The ifS'T images are UV F140PL filter on the left 
and Ha on the right. The aperture used to extract IRS spectra is overplotted on all images. 
An OB association of bright stars are seen on the UV image just outside the aperture box. 
The low- ionization lines are clearly contaminated by these stars. The associated HII region 
is seen on the Ha image. The spectral maps are: [Ne HI] 15.56 /im, [Ne II] 12.81 /im, [S HI] 
18.71 fxm, [ O IV] 25.89 //m [S II I] 33.48 ^m and [Si II] 34.82 //m. The ULX UV counterpart 
identified by lKaaret et al.l (120101 ) has the appearance of a bright star in the ifS'T image, close 
to the center of the MF 16 nebula, and its position is shown with a green X on the spectral 
maps. 
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Fig. 3. — This figure sho ws the SED c o nstru cted in Section 3. The Hubble optical and 
UV data are quoted from iKaaret et al.l (120 lOl ). The base model selected as the best fit 
from the Cloudy simulation (blue) is based on the BMC model fit to the X-ray data and a 
09 I supergiant TLUSTY model star. We also show a model based on the DISKIR model 
(continuous red line). This includes irradiation of the outer disk replacing the companion 
(the disk without irradiation is shown with a dashed red line) The [O IV] edge at 54.93 eV 
is shown as a vertical line. 
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Fig. 4. — Line predictions from CLOUDY simulations. We show the results for both input 
models, BMC and DlSKlR, as model-to-data flux ratios. The perfect match (ratio = 1) is 
shown as a horizontal line. Non-detections are marked with an arrow pointing upward. 
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Fig. 5. — [O IV] line intensity shocks predictions using Mappings III. The shocks plus 
precursor models vary little with the magnetic parameter, and we only show the model 
with a nominal equipartition value (3.23 /iG cm^/^, the dashed line). The horizontal line 
correspond to the brightest pixel in the [O IV] line map in Fig. [2b- The vert i cal lin e marks 
the higher limit of 300 km s~^ for shock velocities measured by iDunne et al.l (J2000l ). 
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Fig. 6. — Ratio line diagnostics comparing the ULXs NGC 6946 X-1 and Holmberg II with 
other extragalactic ionizing sources: AGN, starbursts, etc. Both plots suggest the ULXs are 
very similar to AGNs and LINERs. a) Diagnostic that uses the [ O IV] line nor r aalize d to the 
[Ne II] line. The AGN, starbursts and ULIRG data is taken from lGenzel et al.l (11998). Many 
of the starbursts and ULIRGs are actually upper limits (see Fig. 3 of lGenzel et al.lll998[ ). For 
Holmberg II ULX we show the lower limit because [Ne II] was not detected, b) Diagnostic 
plot that uses the [S III] and [Si II] lines to distinguish between different sources. We plot 
all the data presented by lDale et al.l (120061 ) (their Fig. 5), which includes: Seyferts, LINERs, 
starburst galaxies, H II regions and extranuclear star-forming regions. For Holmberg II 
ULX we show the lower and upper limits because [Ne II] and [S HI], respectively, were not 
detected. 
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Fig. 7. — Line ratio diagnostics from lWeaver et al.l (120101 ) . NGC 6946 X-1 line ratios clearly 
fit within the overlapping high-ionization and low-ionization AGN regime. This is only the 
second ULX to show such line ratios (see also Holmberg II). 



